Efficient generation of graph states for quantum computation 
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We present an entanglement generation scheme which allows arbitrary graph states to be efficiently 
created in a linear quantum register via an auxiliary entangling bus. The dynamics of the entangling 
bus is described by an effective non-interacting fermionic system undergoing mirror-inversion in 
which qubits, encoded as local fermionic modes, become entangled purely by Fermi statistics. We 
discuss a possible implementation using two species of neutral atoms stored in an optical lattice and 
find that the scheme is realistic in its requirements even in the presence of noise. 
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Bipartite entanglement has long been recognized as a 
useful physical resource for tasks such as quantum cryp- 
tography and quantum teleportation. Similarly, multi- 
partite entanglement is an essential ingredient for more 
complex quantum information processing (QIP) tasks, 
and interest in this resource has grown since its con- 
trolled generation was demonstrated in several physical 
systems [l|, |2|- An important class of multipartite en- 
tangled states are graph states. By using vertices in 
a graph to represent qubits, and edges to represent an 
Ising type interaction that has taken place between two 
qubits, the graph formalism gives an effective characteri- 
zation of entanglement by the presence of edges [3 . Spe- 
cial instances of graph states are the resource used in 
multi-party communication protocols, in quantum error 
correcting codes [^ and in one-way quantum comput- 
ing [i|. 

Initial proposals for the generation of graph states in 
physical systems focussed on qubit lattices of fixed ge- 
ometry, where each qubit interacts only with its nearest- 
neighbors [3|. Such a scheme has been experimentally 
implemented in ID with optical lattices of neutral atoms 
via controlled collisions [ij. The graph states generated 
with this method follow the geometry of the lattice, and 
for 2D/3D square lattices they constitute, together with 
single qubit measurements, a universal resource for quan- 
tum computation |5(. However, the generation of more 
complex graphs, where the set of edges does not trans- 
late into a regular arrangement of qubits, e.g. the quan- 
tum Fourier transform graph state, requires the ability to 
pre-engineer a complicated geometry of the qubit inter- 
actions. A simple scheme in which any graph state can 
be created in a set of qubits with a regular fixed geometry 
is therefore highly desirable, and some progress has been 
made towards this with non-deterministic linear optical 
protocols [g. 

In this paper we propose a scheme for efficiently gen- 
erating arbitrary graph states within a linear quantum 
register via an auxiliary entangling bus (EB). The EB is 
a fixed 'always on' system, equivalent, in a specific limit, 
to a non-interacting fermionic system where qubits trans- 
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FIG. 1: (a) Consider a quantum register 7?. which has 3 graph 
qubits in a state |-|-). (b) Two of them are transferred to 
the EB £ where their state is mapped into local fermionic 
modes, (c) £ evolves via H{ for time r, which results in the 
mirror-inversion of the two qubits. (d) The qubits at the 
mirror-inverted location are transferred back to 71, yielding a 
graph state with 2 vertices IG2). Repeating this procedure 
with different qubits allows any 3-qubit graph state to be 
generated in TZ. 



ferred from the register subsequently become encoded as 
local fermionic modes (LFMs) [3 ■ The dynamics of the 
EB over a fixed time then result in the complete mirror- 
inversion of these LFMs yielding robust phases within its 
fermionic state purely as a consequence of Fermi statis- 
tics. Upon transferring the mirror-inverted LFMs back to 
the register the EB has generated entanglement between 
the qubits. For two qubits this procedure performs a 
controUed-CTz (c-CTz) gate, in addition to the inversion J3|, 
as depicted in Fig. Q and this is sufficient to establish an 
edge between the vertices that these qubits represent . 
We examine an implementation of this scheme where the 
EB is an XY spin chain formed from an optical lattice of 
neutral atoms and investigate numerically its fidelity in 
the presence of noise. Finally we describe how, through 
the effective quantum circuit implemented by transfer- 
ring and evolving more than two qubits within the EB 
for the same fixed time, arbitrary graph states of n ver- 
tices can be generated in at most 0{2n) EB cycles. This 
represents an improvement over the 0{n^) steps required 
in a network model composed of two-qubit gates. 

The starting point for our implementation of the EB 
is a ID optical lattice of ultracold bosonic atoms in two 



long-lived internal (hyperfine) states | a) and | b) . The 
dynamics of this system over N sites is then described by 
a two-species Bose-Hubbard model (BHM) given by Q 

n=l ^ ^ 

N~l 

- Y, {t>ian+l + thihn+1 + H.C.) + iJs, (1) 



where an(6„) is the bosonic destruction operator for an 
a(&)-atom in the nth site, and Hb = {B/2)^^{a'lan — 
bl^bn) is the contribution of a uniform external field. The 
parameters t'^^''-' and U°'^''\ U""^ are the laser-intensity- 
dependent hopping matrix elements and on-site interac- 
tions for atoms in states |a)(|6)) respectively. These 
parameters will in general have a spatial profile across 
the lattice. The dynamic controllability and long deco- 
herence times of this system have made it of considerable 
interest for QIP |lfl lllj and for realizing spin models 9] . 
We focus on the two-species BHM in the limit of large 
interactions, U°',U^,U'^^ ^ t°-,t^, which energetically 
prohibit the multiple occupancy of any site. Hopping 
can be then treated perturbatively and to lowest order, 
for an initial Mott insulating state with commensurate 
filling of one atom per lattice site, the effective Hamilto- 
nian is found to be |3] 



Hs- 22 ^n(^nK+l + K 



n=l 



(2) 

Hence we obtain the anisotropic Heisenberg spin model 
in the optical lattice, with \a) = | f) and \h) = \[) 
at each site. The corresponding Pauli operators are 
then al = a^a^ - b^bn, a^ = a'lbn + 6,\a„, crj( = 
— i(ajjfe„ — Hi'^n), while the couplings are given by A^^ = 
(tf + ti^)/i2Uf) - tf/U- - tiyul XI = 4(tf /C/- - 
^flUi) + B/2, and X^J = t^f'jUf. The construction 
of the EB requires the optical lattice parameters to be 
engineered such that U^ = U^ = 2U^^ and i?, = t^, 
thereby ensuring that Af^^ = 0, Af^ = B/2, and that Hs 
reduces to a pure XY spin chain. For simplicity we as- 
sume that only the hopping possesses spatial dependence 
via t'n — ^n — Ty/a^, with a profile obeying < a„ < 1 
so max(t„ ) < r, and that the interaction energies are 
constant over the system as W^ = U^ — 2U'^^ — U with 
both U and T constants. To first order in T/U the dy- 
namics of the optical lattice reduces to an XY spin chain 
with couplings X^ = (2T^/[/)a„. For the moment we 
shall assume ideality where U/T ^ 1; however, later we 
will investigate how the fidelity of this mapping varies 
with U/T and in the presence of noise. 

It is well known that the Jordan Wigner trans- 
formation (JWT) [l3| maps the XY spin chain 
to a non-interacting fermionic Hamiltonian H[ — 



2Xf^y, Un — B, and c„ is a fermionic destruction op- 
erator for the nth site obeying the usual anticommu- 
tation relations. We are particularly interested in the 
angular momentum hopping profile |l3l | given by j„ — 
{J/2)y^n{N — n), which we write as j„ = Wan with 
a„ = 2y/{n/N)[l - (n/N)] and W = AT^/U = JN/A 
so max(j„) < W . In this case the projection of Hf 
onto the single fermion subspace of the lattice, Tii, re- 
sults in a Hamiltonian equivalent to Hi — —JSx + Bl, 
where Sx is the x angular momentum operator for an 
'effective' spin-5 particle, with S = {N — l)/2. The 
single-fermion states { | n) = cj^ | vac) } then correspond 
to the z-angular momentum eigenstates {\S,l)^} of the 
spin-5 particle, with 1 1) = | 5, —S)^ , . . . , | A^) = | 5,5)^. 
The dynamics generated in Tii, when Hi is applied for 
a fixed time t — tt/ J, result in the time-evolution uni- 
tary Ui{t) = cxp (i^s) exp {inSx) composed of an over- 
all phase (pB = —Btt/J for TCi and a rotation of the 
spin-5 particle about the a:-axis by tt. This leads di- 
rectly to perfect state transfer over the lattice [13. The 
action of Ui{t) on the single-particle basis follows from 
its equivalence to the z-angular momentum states where 
exp (iTrS'a;) I 5, Z)^ = exp^inS) \S,—l)^. Thus we find 
that Ui{t) I n) = exp (i(/)i) | n), with the phase (j)i — 7r5 + 
4>B and mirror- conjugate location h — N ~n+l. Choos- 
ing B = SJ is sufficient to ensure that the single-particle 
phase ^1 vanishes. The evolution of the fermionic modes 
cjj then satisfies C/cJ^C/'^ — c\, where U = exp{—iHfT), 
and the dynamics of the system describe the complete 
mirror- inversion of the LFMs. 

Under the JWT the N qubit (or spin) states 
I (71, . . . , gjv), with g„ G {0,1}, of the chain are 
mapped to Fock states of the LFMs as | 91, . . . , gw) = 
{c\)'^^ . . . (c]y)'™ I vac) (using lattice operator ordering) 
which describe the occupancy of the system by quasi- 
fermions. The use of fermionic mode occupancy as a 
basis for quantum computation has been proposed be- 
fore |3] • Complications arise since the accessibility of the 
Fock space of real fermions is restricted by superselection 
rules, and also, in contrast to bosons, it does not pos- 
sess a natural tensor product structure permitting inde- 
pendent operations on each mode, thus making it intrin- 
sically nonlocal [lj|. However, since we have focussed 
on a physical system which m,aps to a non-interacting 
fermionic system, the Fock space is fully accessible, thus 
enabling superpositions of states with different numbers 
of quasi- fermions, which are essential for encoding qubits. 

For systems of identical fermions a bilinear fermionic 
Hamiltonian such as Hi suffices to generate mode entan- 
glement, despite describing a non-interacting system |l5| . 
The entanglement generated by mirror-inversion then fol- 
lows straightforwardly from Fermi statistics through its 
action on Fock states as 
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e |gi,...,9Ar)=e-"^«|gjv,...,gi), (3) 

where Eg is the number of anti-commutations of the op- 
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FIG. 2; For 5 qubits we have (a) the quantum circuit C(5) 
equivalent to the dynamics of H{ for a time r, and (b) the 
fully connected graph state | G5) with 5 vertices generated by 
this circuit if all the qubits are initialized in the | +) state. 
Both the circuit and the resulting graph state generalize in 
an obvious way for more qubits. 



erators cjj required to reestablish a Fock state. Specif- 
Q{Q — l)/2, where Q is the number of 



ically Eq 

fermions, i.e. Q = X]n=i ^Cn = J2n=i 9"' ^^'^ ^o phases 
are only acquired between subspaces with different Q. 
Since Eq. ||2J) is written in terms of Fock states, the 
inverse- JWT removes the antisymmetry when mapping 
back to qubits, while leaving the phases acquired between 
fermion-number (or total magnetization) subspaces un- 
touched. Thus the evolution of the EB over a fixed time 
T is equivalent to a quantum circuit C{N) composed of 
c-CTz gates between all distinct pairs of N qubits followed 
by the inversion operator R, as shown in Fig. E{a) for 
A^ = 5. Usefully, \i N — q qubits are in the state | 0), 
then this circuit reduces to C{q) between the remaining 
q qubits, independent of their locations, followed by the 
full inversion R. 

In passing we note that the above results apply to more 
general settings. Suppose we partition the fermionic lat- 
tice into M equal blocks, each labelled by their cen- 
tral site k, and composed of sites m{k). Within each 
block k we consider an extended fermionic mode fj, = 



E 






defined by a single-particle state (p^ con- 



tained entirely within the block k and symmetrical about 
its center. The dynamics of the EB over some fixed 
time will be equivalent to C{M) R as long as j„ and 
Un are chosen such that the dynamics of Hf performs 
mirror- inversion with respect to fj. . We can equally con- 
sider partitioning a ID continuous fermionic system, de- 
scribed by field operators tp'^ix), and defining extended 
LFMs analogously as fl = J^(j^)dx ip^ {x)tj:'^ {x) , with 
(f)^{x) =0, y X ^ m{k). In this case, harmonic trapping 
V{x) — muP'x^ 12 and Gaussian modes y^ are sufficient 
for mirror- inversion over a time r = tt/w, giving an ar- 
rangement conceptually similar to the ID cold-collision 
proposal in [l£j at the Tonks- Girardeau limit [ij]. 

The generation of arbitrary graph states requires the 
EB to be augmented with a linear register of TV qubits, 
and a transfer process which maps a qubit state into a 
LFM in the EB via (t+ h^ c]j, where (T+ is the Pauli 
ladder operator. The register and the EB are then ini- 



tialized in the states ®^ \ 0)^ and | vac) respectively. The 
scheme begins by choosing a set of register qubits F to 
be the graph vertices, and applying a Hadamard trans- 
formation to each of them: | 0) ^ | +) = (| 0) + | I))/V2, 
as in Fig. ^a) for qubits 1 — 3. A subset S of m of 
these qubits is then transferred to the EB and allowed to 
evolve for a time r, as shown (for tti = 2) in Fig. Q^b) 
and Fig.^c) for qubits I and 2. The LFM qubits within 
the EB at the corresponding mirror-inverted locations S 
are then transferred back to the register, yielding a fully 
connected graph state | Gm) between these m vertices, 
as in Fig. ^d). Such a state is locally equivalent to a 
?7i-qubit GHZ state, as depicted in Fig. EJb) for m = 5. 
Overlap between EB and register graph qubits can be 
avoided by choosing |r| = [iV/2] with locations in the 
first half of the register. 

Our system can generate any graph state of n vertices 
in at most (9(n^) steps by utilizing only the two-qubit 
interaction of the EB to establish each edge individually, 
mimicking the network model of two-qubit gates. How- 
ever, by exploiting the multi-qubit circuit implemented 
by the EB dynamics over the same time r, as in Fig.[5fa), 
our scheme can improve this upper bound. Specifically, 
we proceed iteratively with (; = I by (i) transferring the 
gX\\ graph qubit, and all graph qubits gc > g which will 
connect to g, into the EB; (ii) allowing them to evolve for 
a time r creating a complete set of connections between 
these vertices, c.f. Fig.EJb); (iii) then transferring qubit 
g back to the register while leaving the qubits gc to evolve 
for one cycle longer in the EB, subsequently removing all 
the connections between them; (iv) and finally transfer- 
ring the qubits gc back to the register and repeating step 
(i) with g 1-^ g + 1. Thus, any graph oi g — n vertices 
can be generated in at most 0{2n) uses of the EB. 

In the optical lattice implementation, both the regis- 
ter and the EB are effective spin chains. The register is 
created by turning off the hopping completely |ll| and is 
separated from the EB by a potential barrier (see [13 for 
details of this setup). The transfer process between the 
register and EB can then be accomplished either directly 
via state-independent tunnelling ,I&. .IQj , or through the 
implementation of a swap gate between the register and 
EB qubits jlflj, on a time scale much faster than r. 

We considered some dominant sources of imperfections 
within the optical lattice EB. In particular, we investi- 
gated the fidelity of the two-species BHM to spin-chain 
mapping introduced earlier, for finite U/T. We consid- 
ered a system of size iV = 6 initialized in a product 
state I +) (Xi I 0) (gi I +), and computed the exact time- 
evolution of the two-species BHM using the time-evolving 
block decimation (TEBD) algorithm J20j for varying U/T 
over the appropriate inversion time r. Using the effec- 
tive two-spin density matrix for the end sites, the fidelity 
F was computed with the state | G2) obtained from a 
perfectly implemented XY chain. The simulation results 
in Fig.|2|^a) demonstrate that, as expected, the fidelity in- 
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FIG. 3: (a) The fidelity F of the effective XY spin-chain im- 
plemented by the 2-species BHM with the ratio U/T, for no 
noise A = (■), A = 1% (x) and A = 5% (*). (b) A close-up 
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erated efSciently by using an EB whose dynamics corre- 
spond to a non-interacting fermionic system undergoing 
mirror-inversion. By utilizing an EB which is fixed and 
always on the dynamical control required for QIP tasks 
can be reduced to single qubit operations. Here an im- 
plementation of this scheme using an optical lattice of 
neutral atoms was considered in detail. We also note that 
the EB properties are well suited to solid state systems. 
The fidelity of the optical lattice proposal was examined 
not only for the depth ratio U/T, but also in the presence 
of noise, and found to be both realistic and robust. 
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creases with U/T. Given that r = UNtt/16T^, increas- 
ing the fidelity, by deepening the lattice, comes at the 
cost of longer inversion times. However, Fig. Ofb) shows 
that F > 0.99 even at a moderate ratio of U/T — 26. At 
this depth a A = 514 nm optical lattice of ^^Na atoms has 
r = 9.3 ms, while a A = 826 nm optical lattice of ^^Rb 
atoms has t — 79 ms. These are fast enough for multiple 
EB inversions to occur within the decoherence time of 
the system, which is typically of order of a second |l9| . 

Finally, we investigated the effect of jitter in the lat- 
tice laser intensities. For ^^Rb the laser intensities la 
and lb of the a and b lattices were taken as varying in- 
dependently according to some Wiener noise dWa{b){'t) 
with variance A^: Ia(b){t) = lo + dWa(b)it), where Iq is 
the ideal laser intensity. Such laser fluctuations are then 
non-linearly related to the corresponding fluctuations in 
the hopping t„ ' and on-site interaction U°'^''\U°''' ma- 
trix elements of the 2-species BHM [l9( . We assumed a 
simplified version of this noise in which the fluctuations 
alter the hopping and interaction scalings T and U con- 
tained within the overall scaling W of j„, but not the spa- 
tial profile a„. Despite this restriction, this noise causes 
fluctuations in the inversion time t during the dynamics, 
and also breaks the symmetry required to ensure that 
no u^a^^i or spatially-varying a^ contributions occur. 
In Fig. I^a) the fidelity curves are plotted for A — 1% 
and A = 5% of Iq- For A = 5% the fidelity is seen to 
drop off in deeper lattices due to the cumulative effect of 
noise over longer inversion times. Crucially, the fidelity 
curve suffers only a minimal reduction due to A = 1% 
noise, as in Fig.Ufb), and this represents a realistic value 
for the experimental stabilization of the laser intensity. 

We have shown how arbitrary Graph states can be gen- 



[3; 
[4; 
[s; 

[6; 

ir 

[9; 

[10 

[11 

[12 

[13 

[14 

[15 
[16 
[17 
[18 

[19 
[20 



O. Mandel et al, Nature (London) 425, 937 (2003). 

P. Wahher et al, Nature (London) 429, 158 (2004); 

C.A. Sackett et al, ibid 404, 256 (2000); A. Rauschen- 

beutel et al. Science 288, 2024 (2000). 

M. Hein, J. Eisert and H.-J Briegel, Phys. Rev. A 69, 

062311 (2004). 

D. Gottesman, Ph.D thesis (CalTech, Pasadena, 1997); 

A.M. Steane Phys. Rev. Lett. 77, 793 (1996). 

R. Raussendorf and H.-J Briegel, Phys. Rev. Lett. 86, 

910 (2001); R. Raussendorf and H.-J Briegel, Phys. Rev. 

Lett. 86, 5188 (2001). 

D.E. Browne and T. Rudolph, e- print|quant-ph/0405157| 

S.D. Barrett and P. Kok, e-print quant-ph/0408040 

S.B. Bravyi and A.Y. Kitaev, Ann. Phys. (San 

Diego) 298, 210 (2002). 

M.-H. Yung and S. Bose, Quan. Inf. & Comp. 

4, 17 4 (2004); M.-H. Yung and S. Bose, e-print 

■quant-ph/0407212 

L.-M. Duan, E. Demler and M.D. Lukin, Phys. Rev. 

Lett. 91, 090402 (2003); J.J. Garcia-RipoU and J.l. Cirac, 

New J. Phys. 5, 76.1 (2003). 

J.K. Pachos and PL. Knight, Phys. Rev. Lett. 91, 107902 

(2003). 

D. Jaksch and P. ZoUer, e-print l cond-mat/0410614| 

See for example, S. Sachdev, Quantum Phase Transitions 

(Cambridge Univ. Press, Cambridge, 2001). 

M. Christandl et al, Phys. Rev. Lett. 92, 187902 (2004). 

P. Zanardi, Phys. Rev. A 65, 042101 (2002); L.-A. Wu 

and D.A. Lidar, J. Math. Phys. 43, 4506 (2002). 

V. Vedral, Central Eur. J. Phys. 1, 289 (2003). 

T. Calarco et al., Fortschritte der Physik 48, 945 (2000). 

B. Paredes et al. Nature (London) 429, 277 (2004). 

C. Moura Alves and D. Jaksch, Phys. Rev. Lett. 93, 
110501 (2004). 

M. Greiner et al. Nature (London) 415, 39 (2002); 
W.K. Hensinger et al., ibid 412, 52 (2001); D. Jaksch 
et al. Phys. Rev. Lett. 81, 3108 (1998). 
G. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 



